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A completely positive master equation describing quantum dissipation for a Brownian particle 
is derived starting from microphysical collisions, exploiting a recently introduced approach to sub- 
dynamics of a macrosystem. The obtained equation can be cast into Lindblad form with a single 
generator for each Cartesian direction. Temperature dependent friction and diffusion coefficients for 
both position and momentum are expressed in terms of the collision cross-section. 
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The issue of quantum dissipation, and in particular 
of quantum Brownian motion, is a long standing one, 
which attracts physicists since decades (for general ref- 
erences see and still seems to be unsolved. Its rel- 
evance however is growing, especially in connection with 
decoherence and the relationship between classical and 
quantum description a field which seems to be now 
within reach of experimental tests [Q . The classical un- 
derstanding of the phenomenon is quite well established, 
relying on Langevin or Fokker-Planck equations obtained 
considering a particle typically interacting with a bath 
of independent oscillators, so that most of the research 
work has been influenced by these results, amounting to 
a research for a quantum analog or quantum generaliza- 
tion of these equations. The difficulty lies in the failure 
of a Hamiltonian description for such systems, so that a 
clear quantization prescription is missing and one needs 
a thoroughly quantum mechanical approach. The most 
promising results come from the reduced description of a 
particle interacting with some kind of reservoir, thus im- 
pinging on techniques and results of open quantum sys- 
tem theory. In this respect the property of complete pos- 
itivity (CP) has emerged as a very useful and stringent 
requirement in the study of subdynamics inside quantum 
mechanics . The property of CP amounts to ask that 
the time evolution semigroup lA{t) for the irreversible 
dynamics has the structure ly{{t)g — J2a -^a{t)gK}^{t), 
Cl2a -^a{'t)-^o'{t) — 1), SO that in particular positivity 
is preserved, and it origins from the formal requirement 
that interactionless coupling to an n-level system does 
not affect positivity. Indeed CP appears somehow more 
natural if considered in the context of operations and 
quantum measurement in which it originally appeared in 
physics . According to a famous paper by Lindblad ||] 
under suitable mathematical conditions the property of 
CP allows for the determination of the general struc- 
ture of the generators of irreversible time evolutions, even 
though a thorough understanding of the physical limits 
of validity of this property is still on its way so 
that satisfaction of CP by itself does not ensure physi- 
cally meaningful results, and the connection to realistic 
microphysical models is strongly desirable. In fact re- 
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cent work has stressed the connection between CP of the 
time evolution and weak-coupUng, so that for example an 
uncorrelated statistical operator can consistently be con- 
sidered as initial condition, as well as coarse-graining in 
time [|j . Most of the work carried out at a fundamental 
level regarding dissipative systems follows the Feynman- 
Vernon influence functional formalism (see [p]JlC|] and ref- 
erences therein), first used by Caldeira and Leggett [ pl| , 
who have derived a standard model of quantum Brown- 
ian motion, which however has the drawback that it does 
not ensure the positivity of the statistical operator [ 12|. 
This shortcoming has been cured in various ways [ p^p^ 
and within the independent oscillator model also a posi- 
tive, though not CP time evolution has recently been ob- 
tained . On the side of CP time evolutions very little 
has been done at a fundamental level, thus also leading 
to the conjecture that CP dissipative evolutions could 
not be obtained from first principles , while most ap- 
proaches are either axiomatic or phcnomcnological 
so that few insights about what kind of physical systems 
might consistently be described by Lindblad-type master 
equations can be gained. 

In this Letter we give a derivation of CP quantum 
Brownian evolution based on some recent work on the 
description of coherent and incoherent particle-matter in- 
teraction which has already been successfully ap- 
plied to neutron optics This work relies on the 
introduction of a time scale for the description of the 
subdynamics of slowly varying degrees of freedom in the 
Heisenberg picture, in this specific case the heavy parti- 
cle with respect to the background of thermal particles. 
The interaction is in terms of two-particle collisions de- 
scribed by the full T-matrix and the structure ensuring 
CP naturally arises from the resolvent identity of scat- 
tering theory. 

Derivation and structure of the master equation. Let 
H = Hq + Hni + V he the Hamiltonian of the whole con- 
fined system in second quantization, Hq = J2h ^ha^ah 
describing the particle (either a fermion or a boson in 
the state Uh), matter and V their mutual interac- 
tion. We intend to describe a single particle, so that 
we consider for the total system the statistical operator 
g = J2kh'^k^^'^hQkh, where describes matter, g^h is 
a positive matrix with trace one and we have ahg"^ = 0, 
p™ajj = 0, V/i. In order to consider the subdynamics of 
the microsystem we exploit the following reduction for- 
mula, which connects the expectations of operators in the 
total Fock space Hp with those of operators in the one 
particle Hilbert space 

Tthp (Ag) = ^AhkQkh = Tr^(i) (Ag^ 

h,k 

where A has the typical structure A = ^ a^j^AhkO-k = 
J, aJ^(/i|A|fc)afc. We intend to work on a time scale 
T much longer than microphysical collision time, being 
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only interested in the slow dynamics of the particle, and 
we shall therefore approximate the time derivative with 
the following coarse-grained one: 

= - [gkh(t + r) - gkh{t)\ 

T T 



1 r. 

T 



Tr 



-Hj 



Exploiting the cyclic invariance of the trace we are led to 
consider H-picture operators a\(T)ak{T) , to be evaluated 
on the given time scale using a superoperator formalism, 
so that e.g. TL — j-[H, •], using the integral representation 

ak{T)=e^^ak^ — e'- [z - n)-\k ■ 

J -ioc + q 

Introducing the superoperator T(z) = V + V{z — Ti,)^^V, 
which is the analog of the T-matrix, we have 

(z - ny^ = {z- Hor^ + iz- nor'T{z){z - Ho)'^ 

so that in the considered structure al^{T)ak{T), bilinear 
in the field operators, the emergence of a typically in- 
coherent term having the CP structure KagK]^ bilinear 
in the T-matrix naturally appears, thus confirming re- 
cent phenomenological approaches . Using the fact 
that [H, (I'hO'h] — (the interaction potential being 
bilinear in the particle field operators) the restriction of 
T{z)ak to the case of a single particle may be generally 
written iTiT{z)ak = ^h'^h {'i^z)ah, where (z) is an 
operator in the Fock space of the macrosystem only. This 
matrix, which according to the introduction of the time 
scale should exhibit a slow energy dependence, plays a 
central role, accounting for the peculiarities of the inter- 
action between the particle and the considered medium. 
The master equation we finally obtain has the following 
Lindblad form 



f = -^[H, ^] + i ^ [U,gLi, - i{Ll,U,, g}] , (1) 
where H = Hq + V and V = i(Q -f Q+) with 



{k\Q\h)=Ti-H, [TUEk+ie)g^] 



Ek+Ex-Eh-E^-ie 



being gi™ = '"'4 10 (CI the statistical operator describ- 
ing matter at equilibrium and |A), |^) denoting eigenvec- 
tors of i7,n with eigenvalues E\,E^, while \k),\h) denote 
eigenvectors of Hq with eigenvalues Ek,Eh- A detailed 
derivation of this master equation is given in In 
considering many-particle systems important corrections 
due to statistics of identical particles appear. This will 
not be the case here since we consider a single particle 
distinguished from its environment through the Brown- 
ian limit in which the ratio between the masses is very 
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far from unit. A generalization of the formalism to cope 
with dilute many-particle systems, in which statistic ef- 
fects have been accounted for, has been considered in 
Ref. pl[ , aiming at understanding the transition between 
quantum and classical . 

Completely positive quantum Brownian motion. Fol- 
lowing |lq| we will make the general Ansatz (z) = 
j' (Px J dFy%l}\x)u\{y)t{z,x — y)uh{y)'i}^{x) , where we 
have supposed translation invariance in the interaction 
kernel, and "0^, ip denote field operators for the macrosys- 
tem. Introducing creation and destruction operators 6*, 
6 in the Fock space of the macrosystem, we may write 
'^hi^) = I],,^^]?^fer)V(^)^M- Being interested in local 
dissipation effects we may safely suppose that at least 
far away from the boundaries the system is homogeneous 
so as to use as quantum numbers momentum eigenvalues, 
thus obtaining Tkr,ht,{z) = 5p^+p^^p^+p^i{z,\p^ - p^\), 
depending on the Fourier transform of the interaction 
kernel, and therefore 

Tt{z) = Y.^p„+p,M,b'^i{z. (2) 

We may now insert this expression in (P to evaluate 
the different contributions, starting from the last, typ- 
ically incoherent term. In doing this we consider the 
medium as composed of free gas particles, so that the 
energy eigenstates |A), |^) of may be obtained by the 
repeated action of b\ on the vacuum, and we can write 
|A) = I labeling the different momenta. We there- 

fore simply have (A|&J^5^|A) = 5^.ji.n^ and taking the slow 
energy dependence of t into account the contributions 
with A = ^ cancel out in the master equation and we need 
only consider the primed sum for A 7^ ^. Using, for A ^ 



^'7,m)V "-^Y come to, setting Q^^ = Pf^ Pi; 

1 2e^^„ , i{[p + Q^^]'/2M + ^e,Q,, 



^■.i PP' nP- 2m 2m ^ 2M 2M ^ 



p2 (p Q )2 „,2 (p' + Q )2 ' 
2m 2m ^ 2M 2M 

where M denotes the mass of the Brownian particle, whose position operator is x, while m is the mass of the gas 
particles. The Brownian particle is immersed in a non degenerate gas, so that (rip(l ± n,,)} (n^) (1 ± (n^)) w 
Considering now the quasi-diagonality of the density matrix, linked to its slow variability, we substitute in the T- 
matrix and the denominators p, p' with the symmetric expression ^(p + furthermore we use the variables p^, 
q = Q^^, and put into evidence the ratio a — m/AI between the masses, thus coming to 



47rm X ^ X ^ , 1 

PP 1 



p + p 



1 2M + ie, q 



(1 + a)q + a{p + p')-^ - 2p 

q ^ g 



e^np){p\g\p'){p'\e-^''-'. 



The anticommutator term can be treated in an analogous way, so that the final expression for the dissipative contri- 
butions in (^) becomes, neglecting for simplicity the slow energy dependence of the T-matrix 
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47rm , \t{q)\ 

Q 



pp' 



(l + a)q + a{p + p') ■ ^ - 2p^ ■ ^ 



{l + a)q + 2ap-- - 2p„ • - 



e^nP){p\e\p'){p'\e 



{\P){P19} 



We can now go over to the continuum limit in p^ and q, evaluating the integral with respect to p^ in the case of 
a Boltzmann gas, using (n^) = nA^ exp[— /3(p^/2m)], Am being the thermal wavelength of the gas particles, n their 
density and (3 — l/lk^T) giving the temperature dependence. We thus obtain, in the Brownian limit a 4C 1 



(3h 



iA?_ 



,3 |t(g)|2 j3_„2 



(Fq 



q 



To get to the master equation describing quantum dis- 
sipation we want to extract the temperature dependence 
of this expression, in the limit of small momentum trans- 
fer q. We therefore expand the exponential operators up 
to second order in q, which is also equivalent to keeping 
terms at most bilinear in the operators x and p. Due to 
symmetry properties of the coefficients only terms bilin- 
ear in q and of the form {i denoting Cartesian coor- 
dinates) need be retained, so that we have 



(in 



X 



1 i (3 

-2 [h, [h, q\] + -TTTTT^ [Pf [Pf ^]] + TTTTT t^*' {Pf 



IMP 



h2M 



Let us note that in the derivation important cancella- 
tions and compensations arise between the terms com- 
ing from the anticommutator and incoherent part of (|^), 
necessary in order to obtain the final structure, thus con- 
firming the fact that in quantum theory we cannot have 
separate friction and diffusion terms Decisive for 

the determination of the final structure of the equation 
is also the Brownian limit a = m/M <^ 1. Supposing 
without loss of generality the medium to be isotropic, so 
that q^ = ig^, we can define the following coefficient 



^pp ~ 2 " 



2 7r2m2 



(3n 



d^q\iiq)\'qe' 



(3) 



depending on the collision cross-section through the T- 
matrix and obtain the compact expression 



E 



^[x»,[x„e]]-H^[p,,[p, 



+ ^7[x»,{p„ff}] 



where = {(ih/AMyDpp and 7 = {P/2M)Dpp. 

Exploiting (^) for the T-matrix we simply obtain 
for the potential term in the continuum limit V = 
-"^^ / dMp){p\^(^f{Ep, = 0), n being the density of 
the gas particles, so that it essentially depends on the for- 
ward scattering amplitude f{Ep, 9 — 0) as expected, 
and vanishes if the latter does not depend on energy. The 
complete master equation then becomes 



5 



1= -J[Ho+V,^]-^E[^-[x-^]] 

■i— 1 

- ^ E [p- [p»' ^]] - E t'^- {p^' ^}] • (4) 

i—l i—1 

This is the main result of this Letter, a CP time evolution 
for a quantum Brownian particle derived at a fundamen- 
tal level, using a different, new approach with respect 
to the usual independent oscillator model. The equation 
obtained is translationally invariant and has the correct 
Lindblad form [ p2|j24| . In particular the requirement of 
CP amounts to check that DppDqq > ?i^7^/4, which in 
our case is verified with the equal sign, thus uniquely de- 
termining the different coefficients in a structure like (Q) 
apart from an overall multiplying factor. Let us note that 
the requirement of a stationary thermal equilibrium so- 
lution only determines the ratio between Dpp and 7, and 
also CP simply indicates that the coefficient Dqq should 
be different from zero and within some range, without ac- 
tually fixing it. This explains the wide variety of different 
contributions that have been added to the Caldeira equa- 
tion to make it preserve positive definiteness. The fact 
that DppDqq = ?i^7^/4 has as a consequence the follow- 
ing interesting distinctive feature: in order to write (Q) 
in a manifest Lindblad form, only one generator for each 
Cartesian direction has to be introduced, instead of two. 



In fact using the thermal wavelength \m — /MkT 
associated to the Brownian particle and defining the op- 
erators = ^ (xi + j-^Pi^, satisfying [a 
we can rewrite in the form 



a]] = 



i—l 

+ ^Xhj2[kg^l-i{^V^^,g}]■ (5) 

i—l 

This makes an important qualitative difference with a 
more phenomenological model derived by Diosi |po| , also 
linked to the fact that he obtains an equation with 
the asymmetric expression instead of the momentum 
transfer q. This connection between number of genera- 
tors and relationships among the coefficients in a master- 
equation of the form has not been stressed in the lit- 
erature, even though it provides an important qualitative 
feature, helpful in providing clearcut distinctions. In this 
spirit our work also sheds some light on the recent phe- 
nomenological work of Gao |^ and the subsequent fol- 
lowing debate |^,^ . Gao works from the very beginning 
with a single generator V — fix + ivp and this automat- 
ically leads him to obtain a generalized Caldeira equa- 
tion with Dqq = 7/(8MfcBT), so that DppDqq = h^-f^/i 
is verified. This explains the difference from the diffu- 
sion coefficient Dqq = j /{GMk-sT) found in |jl^. De- 
spite the fact that the coefficients in the master-equation 
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with which Gao starts are actually completely fixed by 
the requirement of thermal equilibrium and his choice of 
a single generator, our work provides some fundamental 
evidence in favor of this structure, giving through (||) the 
quantitative estimate 7 — {(3 /2M)Dnp for the relaxation 
coefficient. The heavily criticized p6| Hamiltonian term 
— ^ {x, p} which Gao obtains, however, does not appear 
in so that no fictitious counterterm is necessary: its 
appearance in rewriting in the form (^) clarifies why 
the initial choice of Gao led to this trouble. 

We have thus obtained a new fundamental derivation 
of CP dissipative evolution, driven by collisions with the 
environment, with temperature dependent friction and 
diffusion coefficients expressed in terms of physical quan- 
tities such as the collision cross-section. The associated 
master-equation has the peculiarity of being expressible 
in Lindblad form with only a single generator for each 
Cartesian direction, thus giving some evidence in favor 
of a recent phenomenological model [Eqj, though being 
deprived of its unphysical features |2^]2^. The underly- 
ing calculations, even though recovering a single particle 
description by tracing over matter, are rooted in a second 
quantization formalism conceived for the description of a 
subset of reduced degrees of freedom slowly varying on a 
given time scale jr7| , p]| . A major extension of this model 
could consist in considering the effect of correlated initial 
conditions on the dynamics and on the property of CP, 
as we intend to do in the future. 

I am very indebted to Prof. L. Lanz who followed the 
whole work and I'd like to thank Prof. A. Barchielli and 
Prof. O. Melsheimer for useful discussions. This work 
was supported by the Alexander von Humboldt-Stiftung. 
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